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Abstract 



Here we consider the 2D free boundary incompressible Euler equation with surface 
tension. We prove that the surface tension does not prevent a finite time splash or splat 
singularity, i.e. that the curve touches itself either in a point or along an arc. To do so, 
the main ingredients of the proof are a transformation to desingularize the curve and a 
priori energy estimates. 
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\0 ! I Introduction 

vq 

Tj" \ In this paper we continue the work in [8] and [9] where we show the formation of singularities 

for the free boundary incompressible Euler equations. Here we prove that in two space 
dimensions the free boundary problem develops finite time "splash" and "splat" singularities 
when surface tension is taken into account (see below, in Section III, the precise definition of 
the splash and splat curves). 
rS ■ In order to describe the evolution of a fluid with a moving domain Q(t) C R 2 , the 2D 

incompressible Euler equations are used: 

(v t + v-Vv)(x,y,t) = -Vp(x,y,t)-(0,l), (x, y) G Q(t) (1.1) 

with the fluid velocity v(x,y,t) G R 2 and the pressure p(x,y,t) G R. The vector —(0,1) 
represents the external gravitational force (the acceleration due to gravity is taken equal to 
one for the sake of simplicity). The free boundary 

dn(t) = {z(a,t) = (zi(a,t),z 2 (a,t)) : a G R} (1.2) 

is smooth and convected by the velocity field 

z t (a,t) ■ z a (a,t) = v(z(a,t),t) ■ z a (a,t), (1.3) 

which is assumed to be incompressible and irrotational 

V-v(x,y,t) = 0, V ± -v(x,y,t) = Q, (x,y)en(t). (1.4) 
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Here we study the relevance of considering the Laplace- Young condition for which the pressure 
on the interface dCl(t) is proportional to its curvature, meaning that the surface tension effect 
is considered: 

~P^ a ^ = 2 \z a (a,t)\3 = 2 K - (L5) 

Above r > is the surface tension coefficient. 

The results in this paper can be shown for three different scenarios: 

1. Q,(t) a compact domain: z(a,t) is a 2-7r-periodic function in a. 

2. Asymptotically flat case: z(a,t) — (a,0) — > as a — > oo. 

3. f2(i) periodic in the horizontal variable: z(a,t) — (a, 0) is a 27r-periodic function in a. 

The problem to study here is the potential formation of singularities for the system (1.1- 
1.5) with smooth interface and smooth velocity field with finite energy as initial data: 

n(o) = n„, dn = {M a ) -a£R}, 

v(x,y,0) = v (x,y), / \v (x, y)\ 2 dxdy < +oo. 

The smooth initial curve zq{o) must satisfy the arc-chord condition: 

\zo(a) -Zo(fi)\ > CAc\a-p\, for all a, ^ € 1, (1.7) 

where cac > is the arc-chord constant. The study of this quantity has been employed by 
other authors to prove local existence (see for example [23], [24]). We will quantify how our 
curve z(a) satisfies the arc-chord condition through the following quantity 

\z{a) — z(a — p)\ 

Throughout the paper we will only focus on scenario 3 for the sake of simplicity. From 
now on, we will denote Qq n [— tt,tt] x R by Oo by abuse of notation (a fundamental domain 
in the period). 

We establish the main result in the paper for the system (1. 1-1.5). 

Theorem 1.1 Consider zq{o) — (a,0) £ H k (T) for k > 5. Then there exist a family of 
initial data satisfying (1.6) and the arc-chord condition (1.7) and a time T s > such that the 
interface z(a,t) G H k (T) from the unique smooth solution of the system (1. 1-1.7) on the time 
interval [0, T s ] touches itself at a single point ("splash" singularity) or along an arc ("splat" 
singularity) at time t = T s . 

These solutions can be extended to the periodic 3-D setting considering scenarios invariant 
under translations in one coordinate direction. In [14], Coutand-Shkoller consider additional 
3-D splash and splat singularities. The case with small initial data was treated by Wu in 
the two dimensional case [25] and the three dimensional case was studied by Wu [26] and 
Germain et al. [16]. 
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For other long time behaviour results see Alvarez-Lannes [3], Castro et al. [10] and the 
references therein. 

In order to prove this theorem we proceed as in [8] and [9]. Using (1.4) it is easy to 
declare that v is harmonic in Q(t). This fact allows us to introduce the moment co(a,t) by 
elementary potential theory as follows: 

v{x,y,t) = — r ——-2 uj{(3,t)d/3, 1.8) 

2vr J R \{x,y) - z{f3,t)\ 2 

where PV denotes principal value at infinity. This moment is also known in the literature 
as the vorticity amplitude. Then the system (1. 1-1.5) is equivalent to the following evolution 
equations which are only written in terms of the free boundary z(a,t) and the amplitude 
u(a, t): 

z t (a, t) = BR(z, uj)(a, t) + c(a, t)z a (a, t), (1.9) 
u) t (a,t) = -2BR t (z,u)(a,t) ■ z a (a,t) - ( — — ) (a,t) + (cw) a (a,t) 



^\d a z\ 

+ 2c(a, t)BR a (z, oj)(a, t) • z a (a, t) — 2(z2)a(a, t) + r 



z aa (a,t) ■ z£(a,t) 
\z a (a,t)\ 3 

(1.10) 



br(z,u) = ^-pv i c ';; yj b " u,((},t)dp, (m 



(for details see for example [12, Section 2]). Above BR(z,uj) is the Birkhoff-Rott integral 
defined by 

]_ py r (z{a,t) - z(P,t))- 
2vr 7 R \z(a,t) - z(/3,t)\' 

and c(a,t) is arbitrary since the boundary is convected by the normal velocity (1.3). 

Local existence in Sobolev spaces was first achieved by Wu [23] assuming initially the 
arc-chord condition. For other variations and results see [15, 21, 7, 27, 24, 11, 20, 13, 22, 28, 
18, 6, 4, 19, 1, 2, 12]. 

The strategy of the proof of the main result is to establish a local existence theorem from 
the initial data that has a splash or a splat singularity (notice that the equations are time 
reversible invariant). Since the curve self-intersects (failure of the arc-chord condition), it 
is not clear if the amplitude of the vorticity remains smooth and the meaning of equations 
(1. 9-1. 10). In order to deal with these obstacles we use a conformal map 

P(w) = (tan(|)) 1/2 , weC, 

whose intention is to keep apart the self-intersecting points taking the branch of the square 
root above passing through those crucial points. Here P(z) will refer to a 2 dimensional 
vector whose components are the real and imaginary parts of P(z\ +iz-2). We also make sure 
that fi(f) U dd(t) do not contain any singular point of the transformation P. Then potential 
theory helps us to get the following analogous evolution equations for the new curve 

z(a, t) = P(z(a, t)) 
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and the new amplitude u: 

zt(a, t) = Q 2 (a, t)BR(z, oo)(a, t) + c(a, t)z a (a, t), 



(1.12) 



'Q 2 {a,t)Cj(a,tf 



ut(a,t) = - 2BRt(z,u)(a,t) ■ z a (a,t) - \BR(z,Cb)\ 2 (Q 2 ) a {a,t) - (- 

^ 4|z a (a, t)\ 

+ 2c(a,t)BR a (z,ui) ■ z a (a,t) + (c(a,t)Q(a,t)) a - 2 (P 2 _1 (5(a, t))) a 

\\z a {a,t)\ j a 

(a,t) ■ z^(a,tY 



+ r(Q 



where 



\z a {a,t)\ 3 



Q 2 (a,t) 



(1.13) 



dP 

dw 



(P'\z(a,t))) 



and HP~ l denotes the Hessian matrix of P^ 1 , which is the i-th (i = {1,2}) component 
of the transformation P^ 1 . 

Here, we choose c(a,t) in such a way that |5 a (a,t)| = A(t). This particular choice of c 
was first introduced by Hou et al. in [17] and was later used by Ambrose [4] and Ambrose- 
Masmoudi [5] . The choice of c implies 



It is easy to check that if we take Q = 1 in (1. 12-1. 13) we recover (1. 9-1. 10). 
We also define the function 

~/ s Q 2 (a,t)oj(a,t) . . .. _ 

(p{a,t) = ^ V ' / \ 1 - c(a,t) z a (a,t) 1.14 
2\z a (a,t)\ 

introduced by Beale et al. for the linear case [7] and by Ambrose-Masmoudi for the nonlinear 
one [5]. This function will be used to prove local existence in Sobolev spaces. 

In the sections below, we show a local existence theorem based on energy estimates. Sec- 
tion III is devoted to provide the appropriate initial data for the splash and splat singularities. 
In Section IV we choose an energy which does not need a precise sign on the Rayleigh- Taylor 
function. In Section V we choose a different energy that involves the sign of the Rayleigh- 
Taylor function and the estimates are uniform with respect to the surface tension coefficient. 
These two energies are based on the ones obtained in the non-tilde domain by Ambrose ([4]) 
and Ambrose-Masmoudi ([6]). 
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The Rayleigh- Taylor function is given by the following formula 



a = ( BRt(z,u) + -^-BR a {z,Co) • z£ + 

-2Vy / ^ I ^rv I 



+ Q 



BR(z,uj) + 



2\z n \ 2 a 



{VQ){z)-zi + {VP; l ){z)-~zi. 



(1.15) 



All solutions that we will consider throughout the paper will have finite energy, as dis- 
cussed in [8]. The system satisfies the conservation of the mechanical energy. We define it 
this way: (not to be confused with the subsequent definitions of some other energies, see 
sections IV and V). 



£s(t) 



v(x,y,t)\ 2 dxdy H — / (^(a, t)) 2 d a zi(a,t)da H — / \d a z(a,t)\da 



a f (t) 



£ k {t)+£ p (t)+£ T (t), 



where z(a,t) = (zi(a, t), Ziipt,, t)), u(a, t) = v(z(a,t),t), and £lf(t) = O(t) n [— it, tc] x E 
is a fundamental domain in the water region in a period, then it follows that the energy is 
conserved. 



d£k(t) 
dt 



n f (t) 



v(x,y,t)(v t (x,y,t) +v(x,y,t) ■ Vv(x,y,t))dxdy 



n f (t) 



v(x,y,t)(-Vp(x,y,t) - (0, l))dxdy 



QAt) 



v(x,y,t)(S7(p(x,y,t) + y))dxdy 



»(«/(*)) 



v(x, y, t) ■ n yds + 



8(0/ (*)) 



v(x, y,t) • n —Kds 



/ \ i \ ^ ±, \ t t f n , no_L/ n 9lz(a, t) ■ d a z (a, t) 

Z2(a,t)u(a,t) ■ o a z (a,t)aa-\ — / u{a,t) • o a z (a,t) j-^ — ^ da 

ZJ-tt \d a z{a,t)\ 6 

(1.16) 

where we have used the incompressibility of the fluid (V-i> = 0) and Laplace- Young's condition 
for the pressure on the interface. Next 



d£ P (t) 
dt 



Z2(a, t)dtZ2(ct, t)d a z\(a, t)da + - / (z2(a,t)) 2 dt( 



z 2 (a, t)d t Z2(a, t)d a zi(a, t)da — / Z2{a,t)d a Z2{a,t)d t zi{a,t)da 

J — 7T 

Z2(a,t)u(a,t) ■ d a z ± (a,t)da. (1-17) 
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d£ T (t) t f w d a z(a,t) ■ d a dtz(a,t) r f w d 2 e z(a,t) ■ dtz(a,t) 

-da = / — — ; -; da 



dt 2J_ n \d a z(a,t)\ 2 7_ 7r \d a z(a,t)\ 



T 



n dlz(a,t) ■ u(a,t) r f n d*z(a,t) • d a z ± (a,t) ± 

da = —- / — — — — u(a, t) ■ OaZ{a, t)da 



\d a z(a,t)\ 2 J_ n \d a z(a,t)\ 3 

(1.18) 

Adding all the derivatives we get the desired result. 

II Properties of the curvature in the tilde domain 

In this section we will rewrite the term corresponding to the curvature K(z(a, t)) in the new 
tilde variables z{a,t). 

We will proceed step by step. Let us recall that the curvature is defined by 

— Zttn ( a ' ' Z a( a i ^) 

\z a (a,t)\ 3 

We begin with the term |z Q (a,t)| 3 . We have that 

\z a {a,t)\ 2 = (d a P{z(a,t)),d Q P{z(a,t))} = (VP(z(a, t)) ■ z a (a, t), VP(z(a, t)) ■ z a (a, t)) 
Since P and P _1 are conformal, by the Cauchy-Riemann equations 

VP(z(a,t)) T VP(z(a,t)) = Q 2 (a,t)Id 2 , 

that implies that 

\z a (a,t)\ 3 = Q 3 (a,t)\z a (a,t)\ 3 

We move to the other term 

{z aa (a,t),z^(a,t)) = (d a (VP-\z{a,t)) ■ z a {a,t)) , (VP~\z{a,t)) • z a (a,t)) x ) 
= (VP _1 (5(a,t)) • z aa (a,t), (X7P-\z{a,t)) • z a {a,t)) L ) 
+ (d a (VP" 1 ^,*))) ■z a (a 1 t),(VP- 1 (i(a,f))-z a (a,()) i ) = W + X 

Again, by the Cauchy-Riemann equations 

W = 1 A z aa (a,t),z a {a,t) L ) 
[a, t) 



Developing the terms in X, we get 



• z a (a, t) — | 



%a(a,t) ■ HP 1 1 (z(a,t)) ■ z a (a,t) 
z£(a,t) ■ ifP 2 _1 (z(a,t)) • z a (a,t) 







where HP i 1 denotes the Hessian of the i-th component of P 1 (i = 1,2). Hence, we can 
write X as 

X = -Za(a,t) ■ i?Pf 1 (5(a,t)) • z a (a, t)VP 2 " 1 (5(a, t)) ■ z(a,t) 
+ z^(a,t) ■ HP^ 1 (z(a,t)) • z a (ot,t)VPi 1 (z(a,t)) ■ z(a,t). 



This means that 

K(a, t) = Q(a, + = + M («> *) 



We will now try to simplify further by exploiting the Cauchy-Riemann equations. We can 
calculate the Hessian and the gradient terms as: 



i M 1 ( 5 ) = *(l^)=*(°) 



Therefore the Hessians are 



»(&) -Q?(6)\ K(6) 



1 " V -9(6) y 2 "V »(&) -9(6) 

Calculating further: 

ziHP^~ Za = + m(&) 2 - 

gHP^Sa = »(&)((4) 2 - {zl) 2 ) - 9(6)(244) 
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X 1 = K(a)K(6)(2(^) 2 ^) + R(a)9(6)((4) 2 *i - (z 2 a ) 2 ~z a ) 
+ 9(a)K(6)(-2zl(zt) 2 ) + ) 2 z 2 - (z 2 ) 2 z 2 ) 

X 2 = H(&)H(&)((£) 2 ig - (5 2 ) 2 z 2 ) + K(a)9(6)(-2^(z 2 ) 2 ) 
+ 9(a)9(6)(2(^) 2 z 2 ) + 9(a)«(&)((4) a £ - (5 2 ) 2 ^) 



This means 



X = Xi - X 2 = ({zif + (zlf){zl (K(a)K(ft) + 9f(a)9f(6)) + 4(»(a)9(6) - 3(a)H(6))) 
= ((^) 2 + ^) 2 )(G(^),5 Q ). 

We can see that 



= »(a)3ft(&)z£ - K(a)^(6)z 2 + S(a)3(&)4 + 9(a)K(6)5 2 
= (G(z),4) 

by the Cauchy-Riemann equations. 

If we take one derivative in space of X, we obtain 



8 a X = {{~z l a f + (z 2 ) 2 )(VG(z) • z a , z a ) + {{ilf + (z 2 a ) 2 )(G(z),z aa ) 
= ((4) 2 + (^) 2 )(VG(5) • z a , z a ) + |S a | 3 £(G(z), z^) 

= ((4)' + (^) 2 )(VG(5) • Za, Z a ) - \z a \ S K^, 

This implies 



K = QK - Q 3 ^- =► K a = {QK) a + ^-(VG(z) ■ z a , ~z a ) - KQ a = (QK) a + M 1 + M 2 

\Z\ | | 

Later, we will see that the Mi is a low order term and can be absorbed by the energy. 

Ill Initial data 

For initial data we are interested in considering a self-intersecting curve in one point. More 
precisely, we will use as initial data splash curves which are defined this way: 

Definition III.l We say that z(a) = {z\ (a), 22 (aO) * s a splash curve if 
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1. z\{a) — a,Z2(cu) are smooth functions and 2ir-periodic. 

2. z(a) satisfies the arc-chord condition at every point except at a\ and cvz, with ol\ < ct2 
where z(a±) = zia^) and \z a {a\)\,\z a {a.-})\ > 0. This means z{a,\) = 2(02), but if we 
remove either a neighborhood of ct\ or a neighborhood of a.2 in parameter space, then 
the arc-chord condition holds. 

3. The curve z{a) separates the complex plane into two regions; a connected water region 
and a vacuum region (not necessarily connected). The water region contains each point 
x+iy for which y is large negative. We choose the parametrization such that the normal 
vector n = ( daZ ^Q°^ z '^^ a ^ points to the vacuum region. We regard the interface to be 
part of the water region. 

4- We can choose a branch of the function P on the water region such that the curve 
z{a) = {z\ (a), Z2(a)) = P{z{a)) satisfies: 

(a) z\(a) and Z2{a) are smooth and In -periodic. 

(b) z is a closed contour. 

(c) z satisfies the arc-chord condition. 

We will choose the branch of the root that produces that 

lim P(x + iy) = - e ~ ln/4 

independently of x. 

5. P(w) is analytic at w and ^{w) ^ if w belongs to the interior of the water region. 
Furthermore, (±7r, 0) and (0, 0) belong to the vacuum region. 

6. z{a) 7^ q l for I = 0, 4, where 

^<°-°>- m^)- 

(nil) 

Moreover, we will define a splat curve as a splash curve but replacing condition (2) by 
the fact that the curve touches itself along an arc, instead of a point. 

Let us note that in order to measure when the transformation P is regular, we need to 
control the distance to the points q l . In order to do so, we introduce the function 

m(q l )(a,t) = \z(a,t) — q\ 

for I = 0, . . . , 4. 

We have performed numerical simulations, as explained in [9] with the following initial 
data on the non-tilde domain: 

2° (a) = a + - (-— - 1.9 ) sin(a) + - sin(2a) + - (- - 1.9 \ sin(3a) 
4 V 2 / 2 4V2 / 
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Zn(a) = — cos(a) cos(2q) -\ cos(3a) 

Note that z (?) = z (— f ) (splash). Instead of prescribing an initial condition for w, we 
prescribed the normal component of the velocity to ensure a more controlled direction of 
the fluid. From that we got the initial u>(a, 0) using the following relations. Let tjj be such 
that V^V = v an d ^(o) its restriction to the interface. The initial normal velocity is then 
prescribed by setting 

u^(a)\z a (a)\ = ^f a (a) = 3 • cos(q) — 3.4 • cos(2a) + cos(3q) + 0.2 cos(4a). 

The reader may easily check that the above z® and z\ yield a splash curve, i.e. the 
conditions in Definition III. 1 are satisfied. See Figure 1. 

0.35 1 i i i i i i i 




-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

x 

Figure 1: Splash singularity. The interface self intersects in a point. 

In order to get an initial data for the splat singularity, one only needs to perturb the 
splash curve so that it z\{ol) = on a neighbourhood of both a = ±5-. The normal velocity 
can be the same since it has the right sign (the one that separates the curve). By continuity, 
the Rayleigh- Taylor function should remain positive. 

For the case where the energy is independent on the surface tension coefficient (see Section 
V), we need the curve to satisfy the Rayleigh- Taylor condition initially. This is always the 
case when the surface tension coefficient is small enough. To illustrate this phenomenon, 
we have plotted in the next figure the Rayleigh- Taylor condition for different values of the 
surface tension coefficient and the initial condition described above. We can see that for 
small enough values of r (0 and 0.1): the Rayleigh- Taylor condition a is strictly positive. For 
bigger values of r, the Rayleigh- Taylor condition a has distinct sign. 

IV Energy without the Rayleigh- Taylor condition 

In this section, we prove local existence in the tilde domain, where the time of existence 
depends on the surface tension coefficient. This theorem has the advantage that the initial 
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Figure 2: Rayleigh- Taylor function for different values of r: r = (blue), r = 0.1 (red), 
r = 0.5 (green), r = 1 (black) 

data does not need to satisfy the Rayleigh- Taylor condition and it works for every r > 0. 

Theorem IV. 1 Let k > 3. Let z° (a) be the image of a splash curve by the map P parametrized 
in such a way that \d a z°(a)\ = -k-, where L is the length of the curve in a fundamen- 
tal period, and such that z^(a),z%{a) G H k+2 (T). Let w(a,0) G H k+ ^(T). Then there 
exist a finite time T > 0, a time-varying curve z(a,t) G C([0, T]; H k+2 ), and a function 
u}(a,t) G C([0, T]; H k+ z ) providing a solution of the water wave equations (1.12 - LIS). 

The proof below is based in the following energy estimates: 



IV. A The energy 

We will define the energy for k > 3 as 



Ei(t) = ££\t)+2\z, 



S |3 



Q 



2k+l 



d k {K)) 2 + l - J Q 2k+2 d k (u)A(d k (u)) + 



2\z a \r 2 



B 



Q 2k+3 (d k a (u>)f 



c 



££ 2 (t) 



l?ll 2 
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where m(q l )(t) = min ae f q l (ct, t) for I = 0, ... ,4 and A = (—A) 1 / 2 . From now on, we will 
denote the Hilbert transform of a function / by H(f), where 



n 2 tan(§ 

Recall that the operator A can also be written as A(/) = d a H(f). 
IV. B The energy estimates 

The energy estimates for EE were proved in [12] and in [8]. In this section we will focus on 
the new terms (A, B and C). 

IV.B.l K 
Proposition IV. 2 

O 2 1 
K t = NICE3 + -^—H(uj aa ) + —^{Q 2 ) a H(u a ), 

Z | Z(x | | Z a I 

where NICE3 means 

j Qid k a {K)d k a {NICE?,) < CE p k (t) 
for some positive constants C,p and any j. 
Proof: We start writing Kt 

Kt = 5 Zfxt ■ z a z aa ■ z^ + I ~ 13 (z aa t ■ z^ + z aa ■ z^t) = -^0 + -Pi + -^2 

I Z<x I I %a I ^ ' 

Calculating further Pq we get that 

P = t^(Q 2 BR + cz a ) a ■ z a z aa ■ ^ = NICE3, 



by the estimates proved in the Appendix. 
On the one hand, developing P2, we obtain 

-^"2 H? ( Zaa ' %at ) — T"^ ( ^aa ' ^at 



,3 l^aa ^at)— I ~ 13 \*aa 



1 /O 2 Z X 

NICE3 - ^--SL-H{u a ) + c a z a )-~zta= NICE3, 
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since c a is as regular clS CO ; "Zqlol 

and therefore bounded in H . On the other, Pi gives rise to 
1 /_ _,\ 1 



-fl = p-f3 ( Zaat • Z a ) = p-fg ((<3 PP)aa + (cZ a ) Qa ) • 5 Q = P^l + Pl >2 
I I ^ ' | I 

We can further develop Pi 5 2 to obtain 

P 12 = NICE3, 

since the terms vanish either by integrating by parts, by being a dot product between two 
orthogonal vectors or because c a = NICE3. We also have that 

P ltl = NICE3 + (2(Q 2 ) a BR a + Q 2 BR aa ) ■ = NICE3 + P 1X1 + P lX2 

I | 

The only term in BR a which is not NICE3 is when we hit with the derivative in Co. Therefore 



P w = NICE3 + -±-2(Q 2 )J-H 



, | & 

Finally, regarding Pi 12 an d keeping in mind that hitting with all the derivatives in z leads us 
to a term which has the factor z aaa ■ z a = — | z aa \ 2 , giving us the extra regularity we needed 
to integrate the term. 

Pi,i, 2 = NICE3 + ^z a ■ J ]s(a) _ m \' u <~WdP) 

= NICE3 + -^—H(u aa ). 

£\z a I 

We should notice that there doesn't appear a term proportional to H(Cu a ) since the kernel 
that results from subtracting the Hilbert transform has room for two derivatives instead of 
one. 

Adding all the previous estimates together we get the desired result. □ 



IV.B.2 Co 

We first notice that M\ (one of the terms in the curvature) is of the order of z a and therefore 
it can be absorbed by the energy. Hence 

d a K = (KQ) a — KQ a + low order terms 

We will follow the proof done by Ambrose in [4]. Taking into account the estimates for 
the implicit operator done in [12], we are left to see the impact of the Q factor in the singular 
term (cco) a , since the impact into the others is either trivial (the ones that come from the 
factor proportional to the curvature) or is zero (the rest of the terms). 
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Lemma IV. 3 

d k (c a u) = NICE35 + 9^H(d k a (K)), 

where NICE35 means 

J Qik{d k a (Cj))NICE35<CEl{t) 
for some positive constants C,p and any j. 

Proof: The most singular term is when we hit all the derivatives in c a , since if we hit all of 
them in oj, that term would belong to NICE35. Developing the new terms 



d k (c a u)) = NICE35 - ud k (Q 2 BR) 



NICE35 - «| . <£ («. ■ / M^fc^ 

Far V J \ z ( a )- z {(3)\ 



NICE35 - ^r~jA (z a ■ H(~z^ a ) 
NICE35 + -^jH (d k (K)\ . 



□ 

Lemma IV. 4 

d k {cQ a ) = NICE35 , 

where NICE35 means 

J Qi 'A(d k {Co))NICE35 < CE p k {t) 
for some positive constants C,p and any j. 

Proof: The most singular term is when we hit all the derivatives in Co a , since if we hit all of 
them in c, that term would belong to NICE35. Thus, we have to estimate 



Q>Hd* + \Z)d* +1 (u,)c 



k+l. 



d k+1 (cj)H(d k+ \Z)Q*c) 



k+l. 



- j d k +\oj) H(d k a +\Cu))cQi -H(d k+1 (Co)cQJ) 



< CE p k (t), 



and therefore it is NICE35. 



□ 
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IV. C Calculations of the time derivative of the energy 

Using the previous lemmas and propositions, we can get the following estimates for the 
derivative of the energy: 

^ = OK + 2 J Q 2k+1 {d k (K^j d k {Q 2 H{u aa ) + 4QQ a H(u a )) 

= OK + 2 J 2kQ 2k+2 Q a (d h a (K))dt 1 (H(£o aa )) 

+ 2 f Q 2k+3 (d k (K))d k (H(Q aa )) 

+ 2 f AQ 2k+2 Q a (d k a {K))d k a (H(u> a )) 
= A 1 + A 2 + A 3 , 

where we will say that a term is OK if it is controlled by the energy. 
We should be careful while estimating Bt because 

^ = OK +i| Q 2k + 2 d k (u }t )A(d k a (Q)) + ~ f Q 2k+2 d k {u>)A{d k (u t )) 
= OK Q 2k+2 d k (oj t )A(d k (u)) + ±- J A(Q 2k+2 d k (ti))d k (u t ) 

= 0K+ lj Q 2k + 2 d k a {u t )K{d k a m + \ j(Q 2k+2 ) a H(d k (u))d k (u t ) 

Hence 

§ = OK + 2 - [ Q 2k+2 A(d k (u))^H(d k (K)) 
at t J 2\z a \ 

+ 2 / Q^A(d k (co))d k ((QK + ^X) a ) 

+ f (2k + 2)Q 2k+1 Q a H(d k (Cj))d k+1 ((KQ)) 

= OK+ 2 -f Q 2k+2 A(d k (u))^-H(d k (K)) 
tJ 2\z a \ 

+ 2 f Q 2k+2 A(d k (co))d k ((QK) a )-2 f Q 2k+2 Q a A(d k a {u))d k {K) 

+ J(2k + 2)Q 2k+2 Q a H(d k (u))d k+1 (K) 
= OK + B 1 + B 2 + B 3 + B 4 
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^ = OK + -i_ f Q 2k +^ 2 d k (oo)d k+1 (K) = OK + c 1 

dt \Za\T J 

IV. D Development of the derivative in B 

We start from the development of B 1 , B 2 , B 3 and B 4 . We trivially have: 

B l = - [ Q 2k + 2 A(d k (Cj))^H(d k (K)) 

T J \Z a \ 

B 3 = -2j Q 2k+2 Q a A(d k (u))d k (K) 
B A = OK - y (2fc + 2)Q 2k+2 Q a H(d k+l (Co))d k (K) 
We now look at B 2 . We can decompose it in the following way 
B 2 = 2 J Q 2k+2 A(d k (u))d k (Q a K + QK a ) 

= OK + 2 J Q 2k+2 A{d k a m{Q a d k a {k) + Qd k a + \K) + kQ a d k a {K^ 
= OK + 5 2 - 1 + B 2 ' 2 + 5 2 - 3 
We can write down the terms B 2,1 and B 2,3 in the form 

B 2 ' 1 = 2 J Q 2k+2 H(d k+ \Co))Q a d k (K) 

B 2 ' 3 = 2kJ Q 2k+2 H(d k+1 (u>))Q a d k (k) 
Integrating by parts in B 2,2 we establish 

-B 2 ' 2 = -2 f Q 2k+3 K{d k+l {Co))d k (K) 

-2(2^ + 3) J Q 2k+2 Q a A(d k (Co))d k (K) 

— ^2,2,1 I ^2.2,2 

Again, 5 2,2,2 can easily be reduced to the canonical form 

£ 2 < 2 > 2 = -2(2k + 3) J Q 2k+2 Q a H(d k+1 (u))d k (K) 

IV. E Collection of the terms 

We will split all the uncontrolled terms into three categories: high order and low order types 
I and II and we will see that the sum of the terms in each category adds up to low enough 
order terms, denoted by OK. 
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IV.E.l High Order 

From A: 

2 f Q 2k ^(d k (K))d k a (H{u: aa )) (A 2 ) 

From B: 

-2 J Q 2k+3 A(d k+1 (u>))d k (K) [B 2 > 2 < 1 ) 

From C: 

No terms from C. 

IV.E.2 Low Order Type I 

From A: 

2 J 2kQ 2k+2 Q a (d k (K)) dt 1 (H(u aa )) (A 1 ) 
2 f AQ 2k+2 Q a (d k (K)) d k (H(u> a )) (A 3 ) 

From B: 

-2 J Q 2k+2 Q a k{d k (u))d k a (K) (£ 3 ) 
- J (2k + 2)Q 2k + 2 Q a H(d k+1 (oj))d k (k) (5 4 ) 
2 J ' Q 2k+2 H(d k+1 (u))Q a d k (K)) (B 2 ' 1 ) 
2k J Q 2k+2 H{d k a +1 {u))Q a d k a (k) (B 2 ^) 
-2(2A: + 3) j Q 2k+2 Q a H(d k+1 (ti))d k (k) (B 2 ^ 2 ) 

From C: 

No terms from C. 

IV.E.3 Low Order Type II 

From A: 

No terms from A. 
From B: 

- [ Q 2k+2 A(d k (u))^H(d k (k)) (B^) 
tJ \z a \ 

From C: 

-J— [ Q 2k+ Wd k (u)d k+1 (k) (c 1 ) 
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IV. F Regularized system 

Now, let z e,<5, ^(a, t) be a solution of the following system (compare with (1.12 - 1.13)): 

z e t ^{a,t) =<t>x*<t> s * {Q 2 (~z e ^)BR(~z £ ^,^)) (a,t) + ^ * (c £ ' 5 ^ * d a ~z £ ^)) (a,t), 

(IV.l) 

2BR t (z £ ' S '",u} £ ' 5 ^) ■ z £ /'" - \BR(z £ ' S ^,oj £ ' S '^)\ 2 (Q 2 (z £S ^)) a 
- (9^!^L) + 27*' S '»BR a (z £ ' 5 'V, oo £ ' 5 ^) ■ z £ />>* + (&> s >i>u?> s >A - 2 ( P^{z £ ^{a, t)) 

+ T f , 5? M '?„ {z^) T HP^z^VP^ ■ - {z^) T HP^~z^VPt ■ 

+T [ Q ^0P~ ) ) " e +» * * ( A( " £ ' V) (IV ' 2) 

z £,s,IJ, (a, 0) = 5o(a) and oj £ ' S ^(a, 0) = (Do (a) for e > 0, <5 > 0, \i > 0. The functions 4>8 and ^ 
are even mollifiers, 



and 

The RHS of the evolution equations for z £,s,tJl and Co £,s ^ are Lipschitz in the spaces i? fc+2 (T) 
and H k+ 2{T) since they are mollified. Therefore we can solve (IV.1-IV.2) for short time, 
thanks to Picard's theorem. 

Now, we can perform energy estimates to get uniform bounds in \i (we just deal with a 
transport term and a dissipative) and we can let \i go to zero. The energy estimates that we 
can get are the following: 



(*) 



| (ll^lfc. + ||^<)||i. + ||^||^ + ± 

< C(S) (\\z^\\%s + ||-F(***")||i- + ll^ll^ + £ ^ii^yj 
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We should note that for the new system without the cf>^ mollifier, the length of the tangent 
vector |S a z' 5 | is now constant in space and depends only on time. Next we will perform energy 
estimates as in the previous case by using the curvature K s from the curve z s . 
Similarly we get (let us omit the superscript 5, e in z s,£ and co s,e ) 
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K t = NICE3 + Jl^fa * fa * H{uj aa ) + -L(Q 2 ) Q< ^ * * H(u a ), 



21? I 3 



d k (c a u) = NICE35 + 9^H(d k (K)), 



d k (coj a ) = NICE35 

and the following collection of terms: 



IV.F.l High Order 

From A: 



From B: 



2 J Q 2k+3 (d k a {K)) d k fa * fa * (H(Co aa )) (A 2 ) 



2 / Q 2fc+3 A(^ +1 (£))<^ * 4>s * d k Q (K) (B 2 ^ 1 ) 



-2-||^ +1 d;||| 2 (D) 



From C: 



No terms from C. 



IV.F.2 Low Order Type I 

From A: 



From B: 



2 J 2kQ 2k+2 Q a (d£(K)) d^fa * fa * (H(Co aa )) (A 1 ) 

2 J AQ 2k+2 Q a (d k (K)) <p s *fa* d k (H(Q a )) (A 3 ) 

-2 J Q 2k+2 Q a A(d k (co))(j)s * fa * d k (K) (B 3 ) 

- J (2k + 2)Q 2k+2 Q a H(d k+1 (ti))fo * fa * d k (K) (B 4 ) 

2 J Q 2k+2 H(d k+1 ("))Q a fa * fa * d k (K)) {B 2 ' 1 ) 

2k J Q 2k+2 H(d k+1 (u))Q a fa * fa * d k (K) (-B 2,3 ) 

-2(2k + 3) J Q 2k+2 Q a H{d k a +l (u))fa * cp 5 * d k (K) (B 2 > 2 > 2 ) 
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From C: 

No terms from C. 

IV.F.3 Low Order Type II 

From A: 

No terms from A. 
From B: 

- ( Q 2k+2 K{d k a (u))^fa*fa*H{d k a {K)) (B 1 ) 

TJ \Z a \ 

From C: 

-J— [ Q 2k+ Wd k (Co)fa*fa*d k+1 (K) (C 1 ) 

We note that throughout this section we have repeatedly used the following commutator 
estimate for convolutions: 

\\h*{d a fg)-g4> 5 *{d a f)\\ L ,<C\\d ol g\\ L oo\\f\\ L ^ (IV.3) 
where the constant C is independent of 5, f and g. 

Also using this commutator estimate we can find all the cancelations we need in the 
previous collection of terms of low order type I and II to obtain a suitable energy estimate. 

Regarding the high order terms, we will do the estimates in detail. We will see the need 
for the dissipative term since there are terms that escape for half of a derivative. 

A 2 + B 2 ' 2,1 + D = 2j Q 2k+3 (d k (K)) fa* fa* (H(d k+2 Co)) 

- 2 J Q 2k+3 H(d k+2 (u))fa * fa * d k (K) - 2e\\d k+1 co\\ 2 L2 

= 2 J d k (K) (Q 2k+ ^ & * fa * H(d k+2 u) -fa* fa* (Q 2k+3 H(d k+2 u))) - 2e 
< ||^J?|| L2 ||a a Q 2fc+3 || LO c||^ +1 w|| i2 -2e\\d k a +l Cj\\ 2 L2 < C(e)EP(t), 

which is uniform in 5. This proves that we can pass to the limit 5 — > 0. 

Finally, by applying the a priori energy estimates to the new system (which only depend 
on e) we can pass to the limit e — > since now we don't have the previous problems and 
A 2 + B 2 ^ 1 = 0. 

V Energy with the Rayleigh- Taylor condition 

In this section, we prove local existence in the tilde domain, where the time of existence does 
not depend on the surface tension coefficient. In this theorem, we need initial data to satisfy 
the Rayleigh- Taylor condition as we explain in Section III. This Rayleigh- Taylor condition 
will hold in particular if the surface tension coefficient is small enough. 
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Theorem V.l Let k > 3. Let z°(a) be the image of a splash curve by the map P parametrized 
in such a way that \d a z°(a)\ = where L is the length of the curve in a fundamen- 
tal period, and such that z\{a),z%(a) G H k+2 (T). Let £(a,0) G H k+ ^(T) be as in (1.14) 
and let u(a,0) G H k ~ 1 (T). Then there exist a finite time T > 0, a time-varying curve 
z(a,t) G C([0,T];H k + 2 ), and functions w(a,t) G C([0, T];H k ~ 1 ) and <p G C([0,T}; H k+ ^) 
providing a solution of the water wave equations (1.12 - 1.13). Assume that initially, the 
Rayleigh- Taylor condition is strictly positive. 

In order to prove this theorem we will use the solutions we have obtained in theorem IV. 1 
for r > 0. We will perform energy estimates on these solutions. 

V.A The energy 

We will define the energy for k > 3 as 



Ei(t) = ss\t) + r^ J Q 2fc+i (d k (K)) 2 + J Q 2k - 2 d k (m(d k a m 



A B 

-\~z a \\ I (C\\K(t)\\ m + K)Q* k ^d k a -\K)k{d k -\K))+2\~z a \ [ C\\K(t)\\ H iQ 2k ~ 2 (d k a {(pf '' 



C D 



m(Q 2k a)(t)' 



E 

where m(Q 2k a) = min ag T Q 2k (z(a, t))a(a, t) and C is a sufficiently large constant such that 
C is strictly positive. Remember that (p was introduced in Equation 1.14. 
At this point is important to notice the following. 

Lemma V.2 The following sentences hold. 

1. Let ip G H 3+ ^, uj G H 2 and z G H k with k > 4. Then u G H 3 . 

2. Let (p G H 3+ ^, uj G H 3 and z G H k with k > 5. Then Q G H 3 5 . 

3. Let Co G H 3+ ^, and z G H k with k > 5. Then ip G H 3 5 . 

This lemma shows that for a fixed r > the energy of this section is equivalent to this one in 
section IV. A. This allows us to use this energy to extend the solutions of the theorem IV. 1 
up to a time T which does not depend on r (for a small enough r). 



V.B The energy estimates 

Again, we will only focus on the new terms (A — E) since the estimates for the other ones 
were proved in [12] and in [8]. 
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V.B.I K 
Proposition V.3 

K t = NICE3B + -^-H(oj aa ) + -J—(Q 2 ) a H(£u a ) 

^ I %a | | Z a I 

= NICE3B + jJ-^H(<p aa ) - rl-{k<p) a , 

I Z a | | Z a | 

where NICE3B means 

' ' Qid k a {K)d k a {NICEW) < CE p k (t) 
for some positive constants C,p and any j. 

Proof: The first equality follows from the proof from the last section since the energies are 
equivalent (see Lemma V.2). We now prove the second one. We begin by using the relation 
(1.14) to get 



We can easily see that 

f — 
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• {Q 2 BR) a = NICE3B 

Za 



since it is at the level of cij a , z aa but we gain one derivative by multiplying by the tangential 
direction. This proves that 

J = NICE3B + ( % 

Fa I \Q 

Looking now to c aa we can see that 

(Q 2 BR) a - • (Q 2 BR) aa = I 1 + 1 2 . 



| Za | | Zql I 

Using the standard estimates, the only thing that causes trouble in l\ is when all the deriva- 
tives hit Cj and therefore 

O 2 

It = NICE3B - K-y— H(u a ). 

1\z a | 

Regarding I2, again, we need all the derivatives to hit BR to get the most singular terms, 
which are 



22 



h = NICE3B 

Q 2 - 



Q 2 



2 - ,~ | 2 



2tt 7„ 7r \z(a) — z 



— — w«(a - p)ap 



(f 



2irJ_, \z(a)-z(fi)\ 2 
1 /' 7r (z aa (a) - z aa (/3)) 



2ir J_ 7T \z(a) — z 



£j aa (a - f3)df3 
Cb a (a - (3)d/3 



NICE3B + 



2Q 2 lza-z- 



| 2 2 



■■-H{u a ) 



Q z a ■ z ai 

12 |5 12 



-H(to a ) 



Q 2 1 u 



\z I 2 2 15 I 2 

I ^ct I ^ I | 



5 Q • H(z aaa ) 



Collecting all the terms from I\ and I2, we obtain 



Q 



- NICE3B + ——H((KCj) a ) 



1 

2\z a \ 
1 



NICE3B +_iJ((^) Q ). 



We can finally write the total contribution as 



■(K<p) a + 



^-cr | 




2(g 2 ) Q 


<■ 


IS I 2 

1 ^"Ct 1 




kQ 2 



Q 2 H ( 4Q a <fc 



Q 3 



Q 2 



<K<p) t 



NICE3B + j—^H (ip aa ) - jz-r(K0) t 



as we wanted to prove. 
V.B.2 <p 



□ 



Throughout this section, we will use the following estimate which was proved in [8] for the 
case without surface tension. The proof is exactly the same for the case with it. 



ip at = NICE2B + 



Q 2 aK + r rp—f ( 

i\Z a 



where NICE2B means 



f QiA(d k a mdt\NICE2B) < CE p k (t) 
for some positive constants C,p and any j. 
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V.C Calculations of the time derivative of the energy 

Using the previous lemmas and propositions, we can get the following estimates for the 
derivative of the energy: 

^ = OK + JL I Q 2k +^d k a {K)d k a (H{y aa )) - r j Q 2k+1 d k (K)d k ((K<p) a ) 

= OK + p- J Q 2k+1 d k (K)d k (H(^ aa )) -t j Q 2k+1 d k a (K)d k a +1 (£)K = OK + A 1 + A 2 

Again, we need to be careful while computing the derivative of B as in Section IV. We obtain 

~ f n 2ifc-2 a lakf,%\\ak-X(,~ \ , f tr&k-2\ utak(~\\ak-l 



2 J Q 2 "- 2 A(d«&))d*-^ at ) + J (Q 2 *-%H(d*(<p))d*-'(v at ) 

= OK - 2 j Q^Aid'mdt 1 
-2 J Q 2k ~ 2 A(d k (<f))d k ~ 1 (Q 2 aK) 
+ r^f Q 2k - 2 A(d k md k (Q 2 (QK) a ) 

- ifr / Q 2k Q a A(d k md k a {R) 

+ / ^(k-l)Q 2k+2 Q a H(d k (v))d k+1 (K) 
= OK + + B 2 + B 3 + 5 4 + B 5 

V.D Development of the derivative of the B term 

We begin noticing that B 1 = OK, as it was proved in [12]. Integrating by parts in B 5 , we 
have that 

B 5 = - [ -^-(k-l)Q 2k+2 Q a H(d k+1 md k (K) 

Furthermore, the only singular terms arising from B 2 are when all derivatives hit either 
K or a, this gives us 

£? 2 = OK -2 J Q 2k A{d k a {(p))d k - l {a)K -2 J Q 2k A(d k (<p))d k ~ 1 (K)a = OK + B 2 ' 1 + B 2 ' 2 . 

However, the only singular term of the Rayleigh- Taylor condition that is not in H k ~ l is 
the one belonging to BRt(z,u>) ■ z a when the time derivative hits u, this means 

B 2 ' 1 = OK - t J Q 2k A{d k {{p))kH{d k {KQ)) 
= - T j Q 2k ^A(d k mKH(d k (K)) 
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Finally, developing B 3 we obtain 

B3 =\h I Q 2k - 2 Md k a md k a (Q 3 k a ) 

\ z a\ J 

+ 1^, / Q 2k ~ 2 A(d k ^))d k (Q 2 Q a k) 

\Za\ J 

Modulo lower order terms we can see that 

B 3 > 2 = OK + -I- [ Q 2k Q a A(d k (^))d k (K) 

We can continue splitting B 3,1 into 

B^OK / g 2fc + 1 A(^( ( ^))^ +1 (^) + ^ / 3kQ 2k Q a A(d k (ip))d k {K) 

I ^tt | J I I J 

= OK --^ / g 2fc + 1 A(^+ 1 (^))^(/?) + ^ l\k-l)Q 2k Q a A(d k (v))d k (k) 
= OK + B 3 ' 1 ' 1 + B 3 ' 1 ' 2 
where in the last equality we have performed an integration by parts. We can observe that 

B 3 ' 2 + B 4 = B 3 ' 1 ' 2 + B 5 = 0, B 3 ' 1 ' 1 + A 1 = 

We will now see that -B 2,2 cancels with the term arising from the derivative of E. Taking 
into account the previous lemmas 



dE 
~dt 



2 J aQ 2k (dtHK)) H(d k+1 (<p)) = OK - B 2 ' 2 



Finally, we will see that the contributions from the time derivatives of C and D cancel 
B 2 ' 1 and A 2 . We start by noticing that, modulo lower order terms A 2 = B 2 ' . Furthermore 

dC f 

— = OK +2r J (C\\K(t)\\ m + k)Q 2k+1 H{d^ + \(p))A{dt\K)) 
^ = OK + 2r J C\\K(t)\\ m Q 2k+1 (d k a {y)) d k+1 (K), 
which, by integration by parts results in 

dC dD i2 21 nrr 
— + — + A 2 + B ' = OK. 
dt dt 

Adding all the contributions, we can bound the derivative in time of the energy by a 
power of the energy. 
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A Helpful estimates for the Birkhoff-Rott operator 

In this Appendix we will prove some of the estimates used throughout the paper for the sake 
of clarity to the reader. 

We begin with a classical decomposition of the Birkhoff-Rott operator. We should notice 
that we can write it in the following ways. On one hand: 



BR(z,uj)(a) = — / , . , -^5 — " in —, - , : u)(a - 3)dB 

V ' ' K ' 2irJ_ n V \z{a) - ~z{a - 2\d a z{a)\ 2 tan(/3/2) J V 

+ — I" ( 10 ^,9 (a) tttttvT ^ w(a - /9)d/9 
27t7„ w V2|3 Q z(a)| 2 tan(/3/2); V P,P 

1 /- f (z(a) - z(a - 0))± \ ^ _ ^ + (q) 



2irJ_ n V \z(a) -z{a-p)\ 2 2\d a z{a)\ 2 tan(/3/2) / v f > r 2\d a z(a)\ 
d^z(a) 



2\d a z(a)\ 2 
On the other hand: 



H(oj)(a) + l.o.t(w). 



w x 1 f n (Ha) -z(a-fi)) L , , „, , . . j „ w(a) P (5(a) - 2(a - /3))" 1 ,„ 



2tt J_ w \z(a) - z(a - 0)\ 2 



2vr \j5(a)-z(a-/3)| 2 4 |z a (a)| 2 sin 2 (f 

w(a) 
2|5 a (a) 



w(a) 



2|z„(a)l 2 



A^a)) +l.o.t(z). 



See [5], [12] for more details concerning the lower order terms. 

We will now prove energy estimates for the Birkhoff-Rott integral, showing that it is as 
regular as d a z. The proof is taken from [12, Section 6]. 

Lemma A.l The following estimate holds 

\\BR{z,u)\\ Hk < CqjF(z)\\ 2 Leo + \\zf Hk+1 + \m 2 Hk y, (A.l) 
for k>2, where C and j are constants independent of z and Co. 
Remark A. 2 Using this estimate for k = 2 we find easily that 

\\d a BR(z,u)\\ L °o < C(||.F(z)||£oo + \\z\\ 2 H3 + M\ 2 H2 y. (A.2) 
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Proof: We shall present the proof for k = 2. Let us write 

BR(z,u)(a,t) = ^f C!(a,P)oj(a - (3)df3 + ^^- H{G})[a) 
where C\ is given by 

U ,P> \z(a) - z(a - P)\ 2 2|d a 5(a)| 2 tan(ft/2)' [ ' 

We shall show that ||Ci||l°° < C\\F{z)\\ 2 Loa \\z\\ 2 c2 . To do so we split G\ = D 1 + D 2 + D 3 
where 

(z(a) - z(a - ft) - d a z{a)l3) L *Lz( n \\ P 1 l 



\z(a) -z(a-f3)\ 2 ' z a ' n \z{a) - z(a - ft)| 2 |5 Q z(a)| 2 ft J 

and 

n r l 1 n 

^3 



|d a 5(a)| 2L ft 2tan(ft/2) J ' 
The inequality 

|5(a) - 5(a - ft) - d Q 5(a)ft| < PIM/?! 2 (A.4) 

yields easily \Di\ < llSlIc 2 !!^)!!^. 
Then we can rewrite D 2 as follows: 

± (d a z(a)(3 - (z(a) - z(a - (3))) ■ (d a z(a)/3 + (z(a) - z(a - ft))) 

|5(a)-5(a-ft)| 2 |cU(a)| 2 ft J ' 

and, in particular, we have 

\d a z(a)P - (5(a) - 2(a - ft))|(|d Q 5(a)ft| + |S(a) - «(a - ft)|) 



l^2| < 



1 5(a) — z(a — (3)\ 2 \d a z(a 



Using (A.4) we find that \D 2 \ < 2||z|| c2 ||J : "(z)|| 2 :oo . 

Next let us observe that since [3 G [— vr,7r] gives \Ds\ < C||J r (5)||i / oo. 
The boundedness of the term C\ in L°° gives us easily 

\\BR(z,u)\\ L2 < C\\T(ml~MhW\v>- ( A - 5 ) 
In d^BR(z, Co), the most singular terms are given by 



1 PV r Q(a 3 M~z{a) ~ d 2 a z{a - ft))- 
-PV .(a - ft) Vz{a) _~ z{a _p )? 



P 3 (a) = -lpV jT £(a - ft) _ _ ^ (5(a) - 5(a - ft)) • (%g( a ) - d 2 a z(a - ft))) rfft. 
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Again we have the expression 

2ir J_ n 2\d a z{a)\ z 

giving us 

|PL(a)| < Cll^ll^pll^dl^lU, + \H(d 2 a Co)(a)\). (A.6) 
Next let us write P2 = Qi + Q2 + Q3 where 



— / (Lu(a — a) — ojya.)) : — r ; — 

2ttJ_^ V H> y " \~z{a) - z(a - /3)| 2 



where A = d a H. 
Using that 

\d 2 a ~z(a) - d 2 a ~z(a - P)\ < \P\ S \\z\\ C 2,5, 
we get |Qi(a)| + |<52(«)| < IMIc 1 ll-^X^IPPII^.i, while for Q 3 we have 



|Q 3 («)| < C\\u\\ l ~\\T(z)\\l-(\\~z\\c* + |A(^)(a)|) 

that is 

\cw 



P 2 (a)\ < (l + |A(^^)(a)|)||cD|| cl ||^)|H|i||^. (A.7) 



Let us now consider P3 = Q4 + Q5 + Qq + Q7 + Qs + Q9, where 

a>(q) (z(a)-z(a-P)-d a z(a)/3) ± . 2 2 m ^ 

Q 5 = — J |z( Q )-z(a-/3)| 4 {{z(a)-z{a-P))-(d a z(a)-d a z(a-P)))dp, 

u(a)d^z(a) r P(z(a)-z(a-P)-d Q z(a)P) ■ (d 2 S(a)-d 2 a ~z(a-(3)) ^ a 
Q6 ~ vr J_ n \z(a)-~z(a-P)\* dP > 

Q7 = - " {a)3 «~ Z{a) d a z(a).r P 2 {d 2 a z{a)-dl~z{a-P)){ } 

vr J_ n "■\z{a)-z{a-P)\ l |<9 Q z(a)| 4 |/3| 4/ 



and 



Qs = - ^ d f[^ d a z(a) ■ (d 2 a z(a) - d 2 a ~z{a - 0))UL - . 1 

vr|<9 a z(a)| 4 J_ 7r v |/3|^ 4sm / (/5/2) / 

LD(a)9^z(a) 2 _ 

g9= -wr aZ(a) ( aZ(a)) - 
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Proceeding as before we get 

|P 3 («)| <C(l + |A(a^)(a)|)||cD|| cl ||^(z)||^||z||^, 
which together with (A. 6) and (A. 7) gives us the estimate 



| (Pi + P 2 + P 3 )(a)| < C(l + \A(d 2 a z)(a)\ + |Pf(^)(a)|)||£|| c i(||j^ 



52 



+ z 



For the rest of the terms in d^BR(z, Q) we obtain analogous estimates allowing us to conclude 
the equality 



□ 



(A.8) 



\\d*BR(z,u>)\\ L2 < C(l + ||^z|| L2 + ||^a)|| L2 )||a;|| 01 11^)11^1 
Finally the Sobolev inequalities yield (A. 9) for k = 2. 
Lemma A. 3 The following estimate will also be helpful 

\\d a BR(z,u) • 8 a z\\ Hk < C(||^(z)||ico + \\zf Hk+ 2 + \\to\\ 2 Hk ) j , 
for k > 2, where C and j are constants independent of z and Cj. 
Proof: In d a BR(z,uj) ■ d a z, the most singular terms are given by 

r> ( \ 1 PT/ f a ~l a , d a ~z(a)-{z(a)-~z{a-fi)) L JQ 

P 2 (a) = ^PV u(a - P ) \~ z{a) _~ z{a _p )? ^, 

P 3 (a) = f Lj( a -P) daZ ^}: [ Z{a l~ Z{a ~f ))L {z{a)-~z{a-P))id a ~z{a)-d a ~z{a-P)))dp. 

vr \z(a) - z(a - P)\ 4 

R2 can be estimated in the same way as P2. Regarding Pi, one can write it as 



Ri(a) = ^PV 



d a z(a) ■ (z(a) — z(a — p)) 1 - d a z a ■ d a z^ 



\z(a) - z{a - p)\ 2 



\z a {a)\' 



dp. 



Now, since d a ui(a — p) = — dpu)(a — (3), one can integrate by parts and bound the resulting 
kernel (which has order -1) giving 

|Pii < c(ii^(5)nicc + ii5|i^ fc+a + \\u\\ 2 Hh y. 

Finally, P3 can be written in the form 



R 3 (a) = PV / uj{a-p) 

7T 



d a z(a) • (z(a) — z(a — p))~ 

\z(a) - z(a - p)\ 4 
d a z{a) ■ d a z(a) 



(z(a) — z(a — p)) ■ (d a z(a) — d a z(a — /?))) 



P\~z a {a)Y- 



d a z(a) ■ d a z(a) 



dp, 



and bound P3 by the kernel (which has order 0) in L°° norm and co in norm. This 
completes the proof. 

□ 
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Then, the following corollary is immediate 



Corollary A. 4 

\\c a \\m < C{\\F{z)\\l°° + \\~zf Hk+ 2 + \M 2 Hk ) j , (A.9) 
for k > 2, where C and j are constants independent of z and Q. 
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